f01qef
f01qef
© Numerical Algorithms Group, 2002.
Purpose
F01QEF Operations with orthogonal matrices, form columns of Q after
factorization by F01QCF or F01QFF
Synopsis
[a,ifail] = f01qef(a,zeta<,ncolq,n,wheret,ifail>)
Description
Q is assumed to be given by
T
Q=(Q Q ...Q ) ,
n n-1 1
Q being given in the form
k
(I 0 )
Q =(0 T )
k ( k)
where
T
T =I-u u ,
k k k
((zeta) )
( k)
u =(z ),
k ( k )
(zeta) is a scalar and z is an (m-k) element vector. z must be
k k k
supplied in the kth column of A in elements A(k+1,k),...,A(m,k)
and (zeta) must be supplied either in A(k,k) or in ZETA(k),
k
depending upon the parameter WHERET.
Parameters
f01qef
Required Input Arguments:
a (:,:) real
zeta (:) real
Optional Input Arguments: <Default>
ncolq integer size(a,2)
n integer size(a,2)
wheret (1) string 's'
ifail integer -1
Output Arguments:
a (:,:) real
ifail integer